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ABSTRACT 

We present the G-virial method which aims to quantify (1) the importance of gravity in molecular clouds in the position-position- 
velocity (PPV) space, and (2) properties of the gas condensations in molecular clouds. Different from previous approaches that 
calculate the virial parameter for different regions, our new method takes gravitational interactions between all the voxels in 3D PPV 
data cubes into account, and generates maps of the importance of gravity. This map can be combined with the original data cube 
to derive relations such as the mass-radius relation. Our method is important for several reasons. First, it offers the the ability to 
quantify the centrally condensed structures in the 3D PPV data cubes, and enables us to compare them in an uniform framework. 
Second, it allows us to understand the importance of gravity at different locations in the data cube, and provides a global picture of 
gravity in clouds. Third, it offers a robust approach to decomposing the data into different regions which are gravitationally coherent. 
To demonstrate the application of our method we identified regions from the Perseus and Ophiuchus molecular clouds, and analyzed 
their properties. We found an increase in the importance of gravity towards the centers of the individual molecular condensations. 
We also quantified the properties of the regions in terms of mass-radius and mass-velocity relations. Through evaluating the virial 
parameters based on the G-virial, we found that all our regions are almost gravitationally bound. Cluster-forming regions appear 
are more centrally condensed. 
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1. Introduction 


Star formation takes place in the dense and shielded parts of 
the interstellar medium. Observations show that the interstellar 
medi um exhibits complicated, irregular, and filamenty y struc¬ 
tures (iGoldsmith et al.ll200^:lMen’shchikov et al.ll201Ql: Schnei- 
der & Elmegreen 1 1979I : I Williams et al.ll20QQh . Theoretically, 
such structures are created by various physical processes such as 
turbulence, gravity, magnetic field, and radiation. 

To understand star formation it is necessary to understand 
how various physical processes affect it. Gravity is a long-range 
interaction, and plays important roles in most astrophysical pro¬ 
cesses at a multiple of physical scale s. This is particula rly true 
for turbule nt molecular clouds fe.g.. iHever et alJl2009l: Kauff- 
mann et al. l2013h . However, characterizing the role of gravity 
on a cloud at various physical scales is not straightforward. 

In observations, the structure of molecular gas can be traced 
by spectral lines in the 3D PPV (position-position-velocity) 
space where molecular gas exhibits complicated structures. 
However, given the rich information obtained observationally, 
few constraints on the role of gravity in the clouds have been 
obtained so far. One major difficulty is to properly quantify 
the irregular structure of the gas. The virial parameter ayt = 
5(T^/?/GM[3 is commonly used to quantify the importance of 
gravity, and to calculate it we must define a region in which cr 


Send offprint requests to: Guang-Xing Li, e-mail: 
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^ Here cr is the velocity dispersion and M is the mass ( Bertoldi & Mc- 
Ke e ll992h. The surface terms and the magnetic terms are neglected. 
See lBallesteros-ParedesI (l2006h for a thorough discussion. 


and M can be evaluated0 As a result, the virial parameter is only 
suitable for the cases where the structures are well defined. The 
morphology of the molecular interstellar medium is generally so 
complicated that in many cases it is difficult to separate individ¬ 
ual objects from a continuous distribution of material. This has 
been further complicated by the fact that we observe the struc¬ 
tures of the sky plane, and the structures identified from the ob- 
serv ations am biased bv vmous projection effects (Beaumont 

et al. l2013l : lDib et al.ll2006l : iPichardo et al.ll200Ql : IShettv et ^ 

[2ml. 

Another difficulty with the virial parameter is that it is lo¬ 
calized. To define the virial parameter we need to define the 
region, and by evaluating the virial parameter of that region, we 
automatically neglect the gravitational interaction between the 
region and its surroundings. As a result, the virial parameter can 
be used to quantify the importance of gravity for the individual 
structures, but it cannot be used to understand the importance of 
gravity for a cloud as a whole. Gravity is a long-range inter¬ 
action, and to understand its importance, we need to understand 
how it works not it one region on one particular physical scale, 
but on a multiple of physical scales. 

One case where such an understanding is required is the 
cluster-forming region. In observations, such regions are found 
to be centrally condensed in both the position-position space 
and the PPV space, and physically we expect that gravity is im¬ 
portant at the centers of the cluster-forming regions and should 
become less important if we move from the centers to the out- 


^ Su ch regions have been de fined with the clump fin d j Williams et al.l 
1 1994) and the dendrogram jRosolowskv et al.ll2008h programs in the 
past. 
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skirts. This stratification, which is of crucial importance for un¬ 
derstanding stellar clustering, awaits to be quantified. 

In this paper we introduce a new method called G-virial 
(available at https://gxli.github.com/G-virial) to quantify the im¬ 
portance of gravity in a variety of situations where mass is traced 
in the 3D PPV space. Instead of dividing the molecular cloud 
into r egions based on iso-in tensity contours in the 3D data cube 
fe.g.. I Williams et aT]ll994l) . the method takes the gravitational 
interactions between all the voxels in the 3D PPV data cube into 
account, and generates maps of the importance of gravity in 3D. 
As a result, it provides constraints on the global importance of 
gravity. This is complementary to methods that calculate the 
self-gravity such as the virial parameter. Such a map of the im¬ 
portance of gravity also enables us to identify gravitationally co¬ 
herent regions in molecular clouds and to quantify the structures 
of the regions. 

In this paper, we introduce the G-virial method (section 
El-Ell, and explore its usage in quantifying the importance of 
gravity and the structures of the molecular condensations seen 
in simulations (section [5]) and observations (section O. We con¬ 
clude in section [71 



2. Problem formulation 

Observationally, molecular gas can be mapped by rotational 
transitions of molecules such as CO. Through proper modeling 
of the emission, it is possible to construct the distribution of the 
molecular gas in the form of a 3D PPV data cube which covers 
a continuous (v, y, v) space 

m(x, y, v) = p(v, y, v) dv dy dv , (1) 

where m represents the mass, v and y represent the spatial di¬ 
mensions, and V represents the velocity dimension. Here p has 
a dimension of M L~^T, and represents the amount of mass per 
unit area per unit velocity. The distribution of p in the v dimen¬ 
sion represents the distribution of gas at different velocities along 
the same line of sight. 

We start with a mass distribution in a 3D PPV data cube (Fig. 
d]). We aim to understand how a particle at voxel i is bound by 
mass from all the voxels j (not excluding /). To achieve this, we 
split our task into two subtasks. First, we need to estimate the 
gravitational boundedness of voxel i by another voxel j in the 
3D PPV data cube. Second, we need to estimate the gravitational 
boundedness of a voxel i based on the information we have about 
all the other voxels j. 

2.1. Boundedness of a voxel pair in a 3D data cube 

We start by defining the gravitational boundedness between a 
pair of voxels (/, j). We are interested in the extent to which a 
particle i is bound to another voxel j. For our purpose, i can 
be treated as a massless test particle, and the boundedness of 
particle i is determined by the mass of the voxel j. 

The available measurable quantities in 3D PPV data cube 
include the mass mj, the spatial separation drij and the veloc¬ 
ity separation dvij. We expect the boundedness to increase with 
nij and decrease with dvij since gravitational attraction increases 
with the mass and decreases with the distance. 

The remaining question is how to make use of the informa¬ 
tion contained in the velocity direction. We expect that voxels 
with large velocity differences are not likely to be bound to each 
other, for two reasons. First, a larger velocity separation im¬ 
plies higher kinetic energy, which consequently decreases the 


Fig. 1: Setup of the problem. We consider a mass distribution 
in a 3D PPV data cube. This is represented as the black grid, 
v-y represents the spatial dimensions and v represents the veloc¬ 
ity dimension. The mass distribution is represented in blue. For 
each voxel /, its G-virial can be determined by taking its inter¬ 
actions with all the voxels j (not excluding i) into account. See 
Sect. Ofor details. 


likelihood of it being gravitationally bound. Second, it has been 
fou nd that the molecular gas follows the Larson’s relation ( Lar¬ 
son 1 198 ih PI where the spatial length and velocity dispersion are 
related by 

6v ~ 6?''^ . ( 2 ) 


As a result, a larger velocity separation implies a larger spatial 
separation, which implies a lower gravitationally boundedness. 

Here, we define the gravitational boundedness contributed 
from voxel j to voxel i as 


= 


SnjSvl ■ 


(3) 


The boundedness increases with mass and decreases with 
spatial and velocity separation. 


Apparent similarities exist between Eq. p and the virial pa¬ 
ramet er. The virial parameter is defined as dBertoldi & McKe j 
Il992lfl 


^^vir ~ c 2 ’ 

5 (t; r 

where m, r, and cr^ are the mass, radius, and velocity disper¬ 
sion of the clumps. The virial parameter is a measurement 
of the gravitational boundedness of an object. The quantity 
defined by Eq. [3] has a dimension of G Mass / Radius 

^ In iLarsonl jl981h the index is 0.38. In iRoman-Duval et~Sl (1201 Ih 
the index is found to be 0.51. 

^ Her e our virial parameter is Fp/2 and in iBertoldi & McK^ 
(Il992h the virial parameter is 2 Fk/Fp. 
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Velocity^, which is the same as that of the virial parameter. 
Therefore, Eq. [3] can be viewed as a generalization of the virial 
parameter to a voxel pair, and a larger Ij^i is related to a larger 
chance for voxel i to be bounded by voxel j. 


2.2. The G-virial 


We define the G-virial of a voxel i as the sum of the gravi¬ 
tational boundedness contributed from all the voxels j (not ex¬ 
cluding i). To be more precise. 


''^G-virial 




Gx 


z 


nijiixi - Xjf + (yi -yj)^) ^ 


= GxJ] 


(Vi - Vj)^ 

mj 


{{xi - Xjf + {yi -yjffivi - Vjf 


(5) 


where 7^^, comes from Eq. |3j and in the second step we take into 
account that the G-virial is measured in the 3D PPV space. 
The physical meaning of Eq. [5] can be understood as follows. 
Eor one voxel /, its G-virial is determined by summing up its 
boundedness with all the other voxels j. 

In Eq. [21 if j is close to i in both the spatial and the velocity 
direction, it contributes more to the G-virial. The contribution 
is proportional to the mass mj, and inversely proportional to dvij 
and Sv^j. 

Eor a continuous distribution of material, Eq. [5] can be writ¬ 
ten as 


t^G-virial(-^9 y? — 


/ 


Gp(x',y',v') 


((x - x')'^ + (y - (v - v')^ 


dx'dy'dv' , (6) 


where the integration is carried out over the whole data cube. 
Here p has a dimension of M L~^ V~^ where M is mass, L is 
size, and V is velocity. 

If the separation between i and j becomes small, the denom¬ 
inator of Eq. [6] approaches zero and the integrand can become 
large. If the density distribution is smooth, the term 


_ p(x,y, v)_ 

((Xi-Xjf + (yi-yjfy^^ ’ 


is still finite and does not induce a singularity. However, the term 


1 

(v - v'Y 

is singular. In our calculations, we change it to 
{v-v’f+cl ’ 

which means we suppress the contribution to the G-virial 
when the velocity separation is smaller than cq, and in our case, 
Co is comparable to the sound speed. We introduce this cutoff in 
order to make the integrand convergent. Nevertheless, there is 
a physical reason behind this. In a molecular cloud, the veloc¬ 
ity dispersion that we observe comes from two parts: one from 
the thermal motion, and the other from the large-scale stream¬ 
ing motion of the gas. The first part is almost scale-independent, 
and the second part increases with the physical size. We are 


mainly interested in the balance between gravity and the large- 
scale streaming motion, and as a result, it is reasonable to discard 
the contributions where 6^ < cq. In our calculations, cq is cho¬ 
sen to be 1 km s“\ which is comparable to the sound speed. 
The effect of changing cq influences the absolute values of the 
G-virial; however, the relative values of the G-virial stay 
unchanged (Appendix [B]). Therefore, the G-virial is a relative 
measure of the gravitational boundedness. In the case where the 
G-virial map of several observations needs to be compared, a 
single value of cq needs to be chosen in advance. 

The reasons for the name G-virial in Eq. [2]include the fol¬ 
lowing: first, we named it G-virial to emphasize its connection 
with the commonly used virial parameter, second, we added the 
letter G to emphasize that our virial parameter is a generalized 
version of the v irial parameter and that it is g lobal. Different 
from t he case of [Bertoldi & McK^ (Il992l) and [Goodman et al] 
([2009[) where the virial parameter is used to quantify self-gravity, 
our G-virial takes all the gravitational interactions between 
gas in the data cube into account. More clarifications concern¬ 
ing the concept of the G-virial can be found in Appendix [Al 


3. Separating components with different veiocities 

One difficulty that we need to deal with is the line-of-sight con¬ 
fusion. When calculating the G-virial, different components 
along the same line of sight are distinguished through their ve¬ 
locities. If two components happen to have the same velocity, 
the G-virial will be over-estimated, since physically unassoci¬ 
ated components will be treated as one single component in the 
calculations. 

We argue that in many cases, the separation in the velocity 
axis provides information so that we can separate different com¬ 
ponents based on the velocity difference or at least alleviate the 
problem. We consider two cases; the first is the line-of-sight 
confusion in a given molecular cloud, and the second is the line- 
of-sight confusion in our Milky Way disk. 

In a molecular cloud, the line-of-sight confusion problem 
can be alleviated if the molecular gas follow s the Larson’s rela¬ 
tion (iLarso nl 198 iLlRoman-Duval et al.l2()TTh because the bound¬ 
edness is proportional to / ~ (Eq. [3]) where and are 

the separations in position and velocity direction, respectively. 
We consider two points along the same line of sight separated 
by ^z- According to the Larson’s relation, the velocity disper¬ 
sion (which is a width of the statistical distribution velocity sep¬ 
aration between the two points) is ~ Therefore, the 

boundedness related to ^z by / ~ ~ 6~^. In other words, the 

contribution is not likely to be large since the velocity is likely 
to be different in an averaged sense. 

In the case of the Milky Way disk, different spiral arms usu¬ 
ally have different velocities, they will be separated easily in our 
calculations. We note that in Eq. [3] the interaction is propor¬ 
tional to 6r~^ and 6v~^. Therefore, a small separation in velocity 
will lead to a much bigger decrease in interaction and finally 
make the contribution to the total G-virial negligible. Since 
the method uses velocity difference to deal with line-of-sight 
confusion, it is more accurate for the cases where the velocity in¬ 
creases with distance (e.g., including systematic expansion and 
contraction where velocity increases with distance). 
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4. Numerical procedure 

In the case of a 3D PPV data cube, for a voxel (x,y,v), the 
G-virial is defined as (see Eqs. [6]and[7]) 


^G-virial(-^? 3^? — 

Gp(x',y\v') 


f 


((x - x')^ + (y- (v - v')^ + Cr 


dx'dy'dv' , (8) 


in the PPV space. Following Equations [3] and [5l we define the 
model G-virial as 


model 

^G-virial, i 


- E 


Gmj 

6rij6vl 


( 12 ) 


where 

Srij = ((X; - Xjf + {ji - yjf + (Zi - Zjf + , (13) 


where the integration is carried out over the whole data cube. 
With an observationally constructed distribution of mass in the 
3D PPV data cube, we can easily construct the G-virial using 
Eq.l 

Equation [8] takes the form of a convolution, and can be con- 
veniently calculated in Fourier space. In real space, the kernel is 


" (v2+y2)l/2(^2+^2) ’ 

and in Fourier space it can be shown tha0 




2 71 

(kl + ^2)1/2 


j:e-‘^o\K\ 

Co 


( 10 ) 


Equation fTOl enables us to calculate the G-virial map ef- 
ficiently. In our calculations, we first make a 3D Fast Fourier 
Transform (FFT) to the observationally constructed mass distri¬ 
bution p(x, y, v) and obtain pk(^x, ky, ky). Then we calculate the 
G-virial in the Fourier space 

ky, ky) = P]^{kx, ky, ky) —Z 2x1/9 ’ (H) 

{kl + k^yi^ Co 

In the last step the G-virial map is obtained by a inverse FFT 
of CKk. 

We found that the absolute values of the G-virial are de¬ 
pendent on Co (see Appendix|B]); however, the relative values of 
the G-virial are insensitive to co. Asa result, the G-virial is 
a relative measurement of gravitational boundedness rather than 
an absolute one. In order for the G-virial maps from different 
regions to be comparable, a unique value of co has to be cho¬ 
sen; Co also needs to be larger than the velocity resolution of the 
observations. 


5. Quantifying the structure of molecular gas with 
G-virial 

We present applications of the G-virial to quantify gravity and 
the structure of molecular gas. To begin with, we apply the 
method to the models where both the 3D density and 3D ve¬ 
locity structure are available. We focus on two aspects: first, we 
present the use of G-virial in quantifying the importance of 
gravity in the 3D PPV data cube; second, we quantify the struc¬ 
tures of molecular condensations bases on the G-virial. 


5.1. Defining G-virial in the RPR space 

To properly evaluate the accuracy of the G-virial method, we 
first define the G-virial”^°^®^ in the PPP space where all the po¬ 
sition and velocity information are available, and then compare 
it with the observationally reconstructed values of the G-virial 

^ Here we follow the non-unitary convention. 


and 


= (Vx - vO" + (v; - + (v^ - vi) 


,J \2 2 

^ ^ model 


(14) 


Here the voxel pair (/, j) has positions ((xi, yi, Zi), (xj, yj, zj)) and 
velocities ((v^, Vy, v^, (Vx, Vy, v^)). The parameter Sx is the spatial 
resolution of the simulation. Writing Eq. [T2]in its integral form, 
we have 


^model ^ f _ Gp{x',y,z') _ 

G-virial, I J _ j,/)2 + _ 2/)2 +^2^1/2 

_ 1 _ 

(Vx(p) - Vx(p')f + (Vy(p) - Vy(p')f 
+ (Vz(p) - Vz(p')f + 

X dx' dy' dz', (15) 


where p stand for (x,y,z) and p' stand for (x',y',z'). Here we 
choose Cmodei = 3 X cq since the denominator of Eq. [15] is com¬ 
posed of a sum of the velocity difference from the three dimen¬ 
sions. Since the G-virial is a relative measure of the gravita¬ 
tional boundedness (AppendixjBj), our results are not affected by 
this. 


5.2. Quantifying the importance of gravity 

We study how the importance of gravity traced in the G-virial 
inferred from the simulated observations is related to the 
G-virial theoretically defined in the models. 

We consider two models. The first (hydrodynamic simula¬ 
tion model) was taken from a numerical simulation of turbulent 
gas. The simulation used here was carried out under periodical 
boundary conditions, and turbulence is inj ected through com¬ 
pressive forcing mode (iFederrath et al.l20()8h . and no self-gravity 
is included. The details of the simulations ar e d escribed in Fed- 
errath et al. (|201Q|) : iRoman-Duval et al.l ([ 2 OIIF FI. Since self¬ 
gravity is not included in the simulation, the line-of-sight confu- 
sion effect is probably over-estimated because, as was found in 
iBeaumont et al.l (l2013h . in the simulations without gravity struc¬ 
tures will overlap more in the PPV space. In this work, we make 
use of the density and velocity cubes from the simulations. 

We make a cutout of a size of 4.8 pc x 4.8 pc x 4.8 pc from 
the snapshot at t = 5 T where T is the crossing timescale of 
the simulation. The computation of Eq. [15] is computationally 
expensive {o{N^)). To save computational time, we choose a 
portion of the simulation with a relatively large gas condensation 
and relatively complicated internal structures (Fig. [2]). 

The second model (the analytical model) is constructed to 
resem ble a typical molecular condensation fe.g. jMcKee & Ta^ 
l2003h . It is defined within Hn < r < rout where r is the radius, 
Tin = 1 pc, and Tout = 10 pc. The density structure is 

P(?') = xpo , (16) 

'•0 

^ The simulations are available at 

http://starformat.obspm.fr/starformat/documentation.jsp. 
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1 PC 


- 2-10 1 2 

X VLSR[kms“^] 

Fig. 2: (a) Map of the peak brightness temperature of the simu¬ 
lated ^^CO(l-O) emission along the velocity axis, (b) Position- 
velocity map of peak brightness temperatures of ^^CO(l-O) 
along the Y-axis (c) Position-velocity map of peak brightness 
temperatures of ^^CO(l-O)) along the X-axis. 

where ro = 1 pc; po = 10^ x where is the mass of the 
H 2 molecule; and jp are chosen to vary from -1.8 to -2.2. The 
velocity structure is parametrized as 

v(r) = X VO , (17) 

ro 

where = 1 pc, is chosen to vary from 0.35 to 0.6, and 
Vo = 1 km s“^. It can be viewed as an expanding sphere. At a 
given radius, the magnitude of velocity is distributed uniformly 
between 0.9 v(r) and 1.1 v(r), and the direction of the velocity 
is distributed uniformly in 3D. The model is sampled through a 
total of 3 X 10^ particles, and then turned into a 3D PPV data 
cube. 

We chose the two models because they represent two typ¬ 
ical situations. In the first case (the hydrodynamic simulation 
model), the model exhibits a high degree of physical complex¬ 
ity. In the second case (the analytical model), the density dis¬ 
tribution is centrally condensed, and the model has well-defined 
mass-radius and velocity-radius structures. We will study the 
behavior of the G-virial in both cases. 

The models are turned into a 3D PPV d ata cube using a sim¬ 
ple ra diative transfer model (Sect. 3.3 of iRoman-Duval et al.l 
I 2 QIII) . It is assumed that the ^^CO( l-O) emission is o ptically 
thin and n l^CQVnfHA ^ 1.7 x 10"^ (iBlake et al.lll987l: Danger 
& Penzias Tl99C)l) . The excitation temperature is assumed to be 
lOK. In reality, the ^^CO(l-O) emission can be optically thick 
and both the abundance and excitation conditions can vary, and 
these can lead to inaccuracies. However, in testing our method, 
we are interested in the cases where mass can be reliably traced 
in the PPV data cube, and observationally this can be achieved 
by modeling several transitions of CO or using more reliable 
tracers. 

In our calculation of the model G-virial 
choose Cmodei = 1 km s"\ and in our calculation of the 
G-virial reconstructed from the simulated map 



0 2 4 6 8 


G-virial Reconstructed 


(a) Hydrodynamic simulation. 



G-virial Reconstructed 

(b) Analytic model. 

Fig. 3: A comparison between the G-virialmodei calculated 
from the model and the G-virialreconstmcted. which is recon¬ 
structed from the simulated observations in the PPV space, (a) 
shows the result for the hydrodynamic simulation and (b) shows 
the result for the analytical model with yp = -2.2 and y^ = 0.5. 
The colors stand for the amount of mass in a given interval. 


we choose cq = 0.3 km s“^ We make this choice because in Eq. 
\T5 \ all three velocity components are considered, and in Eq. [6] 
only one velocity component is considered. Since G-virial is 
a relative measure of the gravitational boundedness (Appendix 
E), the conclusions are not dependent on this. 

In Fig. [3] we compare the G-virialmodei calculated from 
the model with the G-virialreconstmcted reconstructed from the 
simulated observations in the PPV space. The first thing to be 
noticed is that the absolute values of the G-virial differ be¬ 
cause the G-virial is only a relative measure of the gravita¬ 
tional boundedness, rather than an absolute one. 

In general, a larger G-virialmodei is related to a larger 
G-virialreconstmcted- The reconstructed G-virial exhibits a 
higher uncertainty in the hydrodynamic simulation model since 
its flow has complicated structures. 
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Fig. 4: A comparison between between G-virialmodei and G-virialreconstmcted in the PPV space for the hydrodynamical model 
desc ribed in Sect. [5^ In the model, turbulence is injected through compressible forcing, and no self-gravitv is included ( Federrath 
et al. l2008h . The G-virialmnH^i in the PPV space is defined as the mass-weighted average of the G-virialmnH^i in th^PP space, 
which is defined in [121 Here, the grayscale image is the simulated ^^CO(l-O) emission from the model. The red contours stand for 
the G-virialmodeb and the contour levels are (6, 5, 10). The blue contours stand for the G-virialreconstmcted, and the contour levels 
are (3, 6). It worth noting that the G-virialmodei and the G-virialreconstmcted agrees better at regions where mass is traced. This 
can be seen in the v = -0.2 km s“^ channel. They agree better at the regions where mass is traced (e.g., the red arrow), but does not 
agree well at the region without mass (e.g., the blue arrow). 


To look further into the relation between G-virialmodei and 
G-virialreconstmcted, we convert the G-virialmodei from the 
PPP space into the PPV space, and study the connection between 
the two (Fig. n. At a given velocity interval, the G-virialmodei 
in the PPV space is defined as the mass-weighted average of the 
G-virialmodei in the PPP space, and the average is carried over 
the entire line of sight. There are some significant differences 
between the iso- G-virialmodei contours (isosurfaces where the 
G-virial takes constant values) and iso- G-virialreconstmcted 
contours: the iso- G-virialmodei exhibit more structures and 
iso- G-virialreconstmcted is Smoother because the G-virialmodei 
is unaffected by the line-of-sight confusion. Thus, the G-virial 
map constructed in the PPV space is better suited to studying 
the effect of gravity at large scale, but is not suitable for the 
study of the importance of gravity at small scale (e.g., indi¬ 
vidual cores). Although the iso- G-virialmodei and the iso- 
G-virialreconstmcted does not agree completely, they agree bet¬ 
ter at positions where emission is enhanced and mass is traced. 


and tends to differ more where there is little emission. Since we 
are interested in the gravitational boundedness of the gas, this 
disagreement at regions where little mass is present is not as im¬ 
portant as it might appear. 

Since the hydrodynamic simulation model exhibit compli¬ 
cated structures, the line-of-sight confusion has a larger impact 
on the results than the case of the analytical model. However, in 
the gravity-free hydrodynamic simulation mod el the confusion 
probab ly has been over-estimated. According to lBeaumont et al.l 
(l2013h . structures tend to overlap more without self-gravity. 

5.3. Quantifying the internal structures 

For a given region, its G-virial map is centrally concentrated, 
and if we divide the region based on different iso surfaces of 
G-virial values, we get a set of regions that nest inside one 
another. We propose to quantify the structures of the gas con¬ 
densations in the PPV space with the mass-radius relation and 
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velocity dispersion-radius relation derived with the help of the 
isosurfaces of the G-virial. 

Inside a closed region where the G-virial of the voxels is 
larger than a given threshold, we can evaluate parameters such as 
mass, radius, and velocity dispersion. Using a map of mass dis¬ 
tribution m{x, y, v) and a map of the G-virial trc-viriaiC-^, v), 
the corresponding parameters are evaluated by taking all the vox¬ 
els inside a given contour of the trc-viriaiC-^, v) map into ac¬ 
count. The mass of a region is defined as 

m = I p(x,y, v) dxdydv , (18) 

Jr 

where R denotes a coherent region that satisfies 
Q^G-viriai(-^,}^, v) > Qfmin and the Central velocity of a region is 
defined as 


p(x, y, v) V dxdydv 

Vo = —7- • 

p(x, y, v) dxdydv 

The velocity dispersion is defined as 


o-y 


fj^p(x,y, v) (v - vo)^ dxdydv. 1/2 
p(x, y, v) dxdydv ^ 


(19) 


( 20 ) 


The radius of the region is defined by diagonalizing the tensor 
of secon d moments of the posit ion coordinates weighted by the 
intensity (iGoodman et aLluOO^ . and is defined as 


R= 1.91 X (o-minO-max)^^^ , 


( 21 ) 


where (Tminand (Tmax are the dispersions along the major and mi¬ 
nor axes. 

This set of definitions enables us to quantify the structure 
of molecular condensations observed in the PPV space in terms 
of mass-radius and velocity dispersion-radius relations. If the 
molecular gas has a relatively diffuse morphology, we expect to 
see a steeper dependence of the enclosed mass M with the radius 
r; if the molecular gas is centrally condensed, we expect to see a 
shallower M-r dependence. Similar arguments can be applied to 
the dependence of velocity dispersion on radius cTy-r. 

If the molecular gas is centrally condensed, the M-r and cTy- 
r relations of the gas condensations are well defined in the PPP 
(position-position-position) space. It is expected that we can 
reconstruct these relations observationally with the help of the 
G-virial. If a molecular gas condensation has a diffuse and ir¬ 
regular morphology, these G-virial-derived relations can still 
represent the morphology of the gas condensations. The differ¬ 
ence between different regions can be represented in the M-r and 
(Ty-r planes. 

To demonstrate the diagnostic power of the M-r and CTy-r re¬ 
lations, we generated a set of models with different density and 
velocity structures, and compare them with the M-r and CTy-r 
relations reconstructed with the G-virial (Fig. [5l). The dif¬ 
ferences in structure among different regions can be accurately 
reconstructed with the help of the G-virial. Furthermore, the 
slope of the M-r and CTy-r relations are accurately reconstructed 
with the G-virialQ. 

^ To make the G-virial-derived relations overlap with the model val¬ 
ues, we scaled the estimated radius by a factor of 1.4 and the estimated 
velocity dispersion by a factor of 1.45. This adjustment is made in or¬ 
der to make it easier to compare different curves. In practice, the radii 
of the estimated regions are based on the tensor of second moments of 
the position coordinates, and the radii do not necessarily coincide with 
the radii defined in the model. A similar argument also applies to the 
estimation of the velocity dispersion. 




Fig. 6: G-virial as a function of radii for models with different 
yp and yy. 


We also plot the dependence of G-virial as a function of ra¬ 
dius for various models (Fig. Ell. For all the models, we found an 
increase of G-virial towards the centers of the regions. Larger 
G-virial are found for the models with more condensed struc¬ 
ture and lower velocity dispersion. 

5.4. Short summary 

In this section we demonstrated the usage of G-virial to quan¬ 
tify the structure of molecular gas condensations. We summarize 
our findings and remind the reader of the caveats. 

We found that the G-virial can be used to quantify the im¬ 
portance of gravity and to quantify the structures of the gas con¬ 
densations. 

1. The G-virial measured from the PPV data cube is pos¬ 
itively related to the G-virial found in the model (PPP 
space), and a larger G-virial is related to a larger proba¬ 
bility of being gravitationally bound by the ambient gas. 

2. The G-virial map can be used to derive mass-radius and 
velocity dispersion-radius relations. The G-virial-derived 
relations are good reconstructions of the structures of the 
models when the models are symmetric, and can also pro¬ 
vide good representations of the structures if the models are 
irregular. 
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Fig. 5: Comparison between the structure of the model and the structure reconstructed through the G-virial for different sets of 
model parameters. The upper panels show the results where the analytical model is constructed with different 7 p, and the lower 
panels show the results where the analytical model is constructed with different yy. The left panels are the mass-radius relations and 
the right panels are the velocity dispersion-radius relations. Here the radius is scaled up by a factor of 1.4 and the velocity dispersion 
is scaled up by a factor of 1.45 for an easy comparison. 


Several caveats need to be noted. First, the G-virial is 
a relative measure of the gravitational boundedness. A larger 
G-virial means a larger change of being gravitationally bound. 
The absolute values of the G-virial cannot be directly used to 
tell if a region is gravitationally bound or not. The absolute val¬ 
ues of the G-virial are dependent on cq. Therefore, in order for 
the values of G-virial from several maps to be comparable, a 
unique value of cq needs to be chosen in advance. This value 
should be comparable to the sound speed of the medium, and 
in this work we have chosen it to be cq = 1km s“^ We also 
note that the G-virial is reconstructed from the PPV space 
where the projection effect tends to smear out structures; it is 
more suitable for the study of the gravitational boundedness of 
the structures larger than the individual gas condensations. 

Second, the G-virial method can be used to derive the 
mass-radius and velocity dispersion-radius relations for molecu¬ 
lar condensations. However, the results need to be properly inter¬ 
preted. If the structures of the gas condensations are regular and 
are close to being spherically symmetric, the G-virial-derived 
mass-radius and velocity dispersion-radius relations should be 
good reconstructions of the real 3D structure. If the structure of 


the gas condensations are irregular, the G-virial can be used 
to derive these relations. However, in this case, the G-virial- 
derived relations should be viewed as representations of the com¬ 
pactness of the structures in the PPV space, and are not neces¬ 
sarily related the compactness of the structures in the PPP space. 

Finally, we would like to emphasize that the G-virial 
methods works in the PPV space, and like other methods, its 
accuracy is limited by line-of-sight contamination. The method 
relies on using the velocity difference to alleviate the line-of- 
sight confusion. This requires different gas components to have 
different velocities, which is satisfied in many cases. For a sin¬ 
gle object, this requires the velocity separation to increase with 
radius. Whereas this generally holds for individual molecular 
clouds and clumps, it does not hold for some cores or disks 
where the velocity decreases with radius (e.g., collapsing cores). 

6. Applications to moiecuiar ciouds 

In this section we present applications of our method to several 
molecular clouds. Our method requires a map of the distribution 
of molecular gas in the 3D PPV space. Observationally, this can 
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be conveniently achieved using rotational transitions of the CO 
molecule. 

In the simplest case, this can be achieved by observing the 
^^CO(l-O) transition alone. Assuming an excitation temperature 
of 10 K, a abundance of X(i2C0/H2)=8 x 10“^ and a 

to ratio of R(^^CQ/^^CQ)=45, the column density c an be 
estimated as (iRoman-Duval et al]l201Ql : [^mon et al.ll200ll) 

NiH2) = 4.92 X 10^° X — X — L- , (22) 

K km s“^ 

where Tnib is measured in K and Sy is the velocity resolution of 
the data cube measured in km s“^. Changing the excitation tem¬ 
perature to 20 or 3 0 K decreases the derived mass by 40% and 
92%, respectively (ISimon et al.ll200ll) . For an observed ^^CO(l- 
0) data cube rnib(x,y, v), the mass distribution m(x,y, v) can be 
obtained as 


m{x, y, v) 
8 


= 4.92 X 10 


20 


-1 


K 


km s 


mu. 


(-) 

'^cm^ 


(23) 


where is the voxel size of the spatial dimension in cm and Sy 
is the channel width in km s“^. The values of m and p in Eq. [6] 
are related by m = p dx dy dv. 

We apply our G-virial meth od to the publicly available 
data from the COMPLETE survey (iRidge et al.ll200^ . The ob¬ 
servations have a spatial resolution of 46" and a velocity resolu¬ 
tion of 0.067 km s“^ the mean RMS per channel is ~ 0.33 K in 
terms of and the beam efficiency is ~ 0.5. 

With a distance of 250+50 pc and a total mass of lO^M© 
(lEnoch et al .1120061) . the Perseus molecular cloud is among the 
best-studied molecular clouds in the Milky Way. It is composed 
of several distinct regions: Bl, B3, NGC1333, and IC348. The 
Ophiuchus mol ecular cloud has a d istance of 125pc and a mass 
of 7 X 10^ Mq (Ide Geus et al.lll98^ and contains clouds L1688 
and L1689. In particular, L1688 and NGC1333 are cluster¬ 
bearing regions. IC348 also hosts a star cluster; however, it 
is generally considered older than the other regions in Perseus 
dGutermuth et al.|[200^ . 

In this section, we provide a study of the structure of the 
Perseus and Ophiuchus molecular clouds with our G-virial 
method. We present maps of the G-virial, together with an 
analysis of the structure of the regions in the clouds based on the 
G-virial method. 


6.1. Maps of the G-virial 

In Fig. [7] we channel maps of the G-virial in the 3D PPV 
space. Three-dimensional renderings of the maps are presented 
in Appendix O For both clouds, the ^^CO(l-O) maps contain 
complicated and filamentary structure, whereas the G-virial 
maps are smooth and contain fewer structures. Interestingly, 
higher values of G-virial are reached only at cluster-bearing 
regions such as NGC1333, IC348, and L1688, and this high¬ 
lights the importance of gravity in such regions. Different re¬ 
gions are indicated in Fig. [S] 

6.2. Identification of regions 

The G-virial offers a new way to divide the molecular cloud 
into regions. Previously, this has been done either by visual in¬ 
spection, or with contour-based algorithms such as clump find 




RA. (Degree) 


Fig. 8: Regions identified by applying the dendrogram al¬ 
gorithm to our G-virial map. The results from the Perseus 
molecular cloud (upper panel) and Ophiuchus molecular cloud 
(lower panel) are presented. The grayscale images correspond 
to the velocity-integrated ^^CO(l-O) emission and the contours 
correspond to the projected boundaries of the identified regions. 
These regions correspond to the “leaves” of the dendrogram. 
The conventional names of the regions are labeled. 


(I Williams et al.l Il994l) . It is worth mentioning that other 
region-finding algorithm s are available , such as den drogram 
(Ififfisolowskv et al.ll2008h and do champ (IWhitingll2012l) . 

Many of the algorithms (e.g., clumpf ind and dendrogram) 
are contour-based, and they tend to assign voxels to regions 
based on iso-column-density contours in the 2D case or on iso¬ 
intensity contours in the case of a 3D PPV data cube. Molecu¬ 
lar clouds are characterized by a set of complicated hierarchical 
structures, and a naive application of the clumpfind tends to 
produce hundreds of clumps for one single molecular cloud. It 
has also been pointed out recently that contour-based methods 
suffer from superposition and conf usion when the volume filling 
factor of emitting material is large (iBeaumont et al.ll2QT^ . 

In our analysis we are interested in finding gravitation¬ 
ally coherent regions. This can be achieved by applying the 
dendrogram to the G-virial map, which is a measurement of 
the gravitational boundedness. Figure [8] shows the regions iden¬ 
tified by applying the dendrogram algorithm to our G-virial 
map. The identified regions correspond to the “leaves” of the 
dendrogram. In the dendrogram calculations, the minimum 
difference between different contours is set to 0.4, and the re¬ 
gions correspond to the leaves of the dendrogram. Each identi¬ 
fied region consists of a continuous list of voxels in the 3D PPV 
space. In Fig. [8]the projected boundaries of the regions are plot- 
tedS. 

^ All the regions were originally defined in 3D. 
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Fig. 7: Channel maps of the ^^CO(l-O) emission and the corresponding G-virial. The upper panel shows the results from the 
Perseus molecular cloud and the lower panel shows the results from the Ophiuchus molecular cloud. Grayscale images stand for 
^^CO(l-O) emission and the red contours represents the G-virial. Contour levels start from 1.2 in steps of 0.8. 
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Compared to clumpfind, a combination of the G-virial 
method and the Dendrogram method tends to identify regions 
that are coherent while the clumpfind tends to break up those 
regions. The reason is that our G-virial maps are much 
smoother than the original map0. Since the G-virial maps 
are much smoother, the results obtained with the G-virial 
maps are much less sensitive to the technical parameters (e.g., 
the minimum separation of the contours) used in the computa¬ 
tions compared to the case with the intensity maps. Therefore, 
the G-virial offers a more robust definition of regions. With 
G-virial we can study the structure of molecular gas in terms 
of these gravitationally coherent regions. In the 2D case, our 
method is similar to the gridco re code which is based on 2D 
projected gravitational potential (iGong & Ostriked[201lh . 

6.3. Internal structure of the regions 

In this section we study the internal structure of the molecular 
gas in the individual regions. Previously, the structure of molec¬ 
ular gas was studied in terms of the clumps, and quantified using 
parameters that are evaluated for the whole clumps. While these 
methods do provide constraints on the role of gravity inside the 
regions, the role of gravity on scales larger than the individual 
gas condensations is neglected. The spatial structure of molec¬ 
ular gas larger than the clump scale can be quantified using the 
dendrogram algorithm. However, the dendrogram method is 
contour-based, and only the connection between the adjacent re¬ 
gions in the tree diagram are preserved. This makes it difficult 
to use dendrogram to quantify the structure of molecular gas on 
the large scale. 

Here we provide an analysis of the structure of the individual 
regions with the G-virial method. As we have shown in sec¬ 
tion 15.31 the intrinsic mass-radius and velocity-radius relations 
in the model can be reconstructed with our G-virial method 
to a good accuracy if the velocity dispersion increases with the 
size. 

In Fig. [9] we plot the dependence of velocity vq and veloc¬ 
ity dispersion cr^ on the radii of the regions. The inner parts 
of the regions are clearly at the center of the outer part of the 
regions in velocity space. This means the central part of the re¬ 
gions remains quiescent with res pect to the outer pa r ts, and this 
is con si stent with the findi ngs bv l Andre et al.l (l2007l) : [Kirk et al.l 
(l2010h : IWalsh et alJ (12004 . 

We assess the importance of gravity in these regions with the 
G-virial method. In our method, a larger G-virial is related 
to a larger importance of gravity, and as a result the importance 
of gravity at a given region can be quantified by measuring the 
amount of gas at different thresholds of the G-virial parameter. 
In Fig. [TOl we plot the dependence of gas mass on the G-virial 
threshold for different regions, and in Fig. [TT]we plot the depen¬ 
dence of the G-virial on the radius. For all the regions, the 
G-virial increases towards the centers. For cluster-bearing re¬ 
gions such as NGC1333 in Perseus and LI688 in Ophiuchus, 
much higher G-virial values have been reached at their cen¬ 
ters. This implies that gravity is more important for the clustered 
mode of star formation. 

With our method it is straightforward to derive relations such 
as the mass-radius and velocity-size relation. Recently, there 
has been growing interest in quantif ying cloud structure in ter ms 
of various mass-size relations fe.g.. lKauffmann et al.ll2010j^ . 

^ This is similar to the case of ISmith et all j2009l) where they identified 
structures form a simulation in the 3D position-position-position space 
based on gravitational potential. 





H NGC1333 

]—[ IC348-B3 
H L1688 
H L1689 


10 " 

Radius [pc] 


Fig. 9: Dependence of velocity and velocity dispersion as a func¬ 
tion of radius. Here the solid lines represent the velocities of the 
regions, and the bars represent the velocity dispersions of the re¬ 
gions at given radii. The vertical extent of the errorbars are the 
3D velocity dispersion = ^^3 cr^. 



Fig. 10: The amount of molecular gas enclosed in regions with 
different thresholds of the G-virial. The x-axis is the threshold 
of G-virial which we take to define the region, and the y-axis 
is the amount of gas that is enclosed in the region. 



Fig. 11: G-virial as a function of radius for different regions. 
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Fig. 12: Enclosed mass as a fun ction of radius for different re- 
gion s. The scaling relation fromiLarsonI (1 19811) and Kauffmann 
et al. (l2010bl) are added. For all the curves, the G-virial de¬ 
creases with increasing radii. 


Fig. 13: 3D velocity dispersion as a function of r adius for differ- 
ent regions. The velocity-linewidth relation from iLarsonI (Il98lh 
is also plotted. 


In Fig. [12] we plot the mass-size relation obtained using our 
G-virial method (see Section[5]). Our mass- size relation is dif¬ 
ferent from that of I Kauffmann et af] (l2010bl) since we use the 
G-virial contours to define our regions, and mass is evalu¬ 
ated within a giv en region where the G-vi rial is larger than 
the threshold. In iKauffmann et al.l (1201 Obi) the mass in evalu¬ 
ated within a region where the column density is larger than a 
given threshold. In our case, along one single line of sight, only 
the gas that stays within a given surface of a constant value of 
the G-y irial is taken into account while in IKauffmann et all 
(l2010bl) al l the mass along the line o f sight is taken into account. 
Similar to IKauffmann et al.l (l2010bl) . we also found that cluster¬ 
bearing regions such as LI688 and NGC1333 are more massive 
in terms of molecular gas than the regions devoid of clusters at a 
gi ve radius. In Fi g. [12] we also plot the scaling r elations proposed 


m 


Larsonl(ll98ll. m{r) = 460Mo(r/pc)^-^), and IKauffmann et aP 


(l2010bL m{r) = 400Mo(r/pc)^-^) , and found that both provide 
approximate descriptions of the structure of the regions. 

Molecular clouds are believed to be turbulence-dominated, 
and this is mainly inferred fro m the fact tha t the clouds obey 
the velocity-linewidth relation dLarso l^ ll98lh . In Fig. [B] we 
plot the 3D velocity dispersion of our regions as a function of 


radius. Here the 3D velocity dispersion is V3 times the velocity 
disper sion evaluated in Eq. [20] The scaling relation from lLars^ 
(Il98ll) is also plotted (L/pc ~ 1.01 x (cTy/ km s“^)^-^^). In 
our case, the velocity-linewidth relation is evaluated for centrally 
condensed objects, and a smaller physical scal e is rela t ed to t he 
inner part of a gas condensation, whereas in iLarso BdMIl) a 
smaller phwcal scale is related to a smaller size of a subregion 
in a cloud E . For our cent rally condensed objects, the scaling 
relation from lLars^ (Il98ll) seems to be valid. We also note that 
cluster-bearing regions such as NGC1333 and LI688 have much 
larger velocity dispersi ons at large radii. 

Heyer et al. (l2009l) found a dependence of the scaling coef¬ 
ficient (Ty/r^/^ of Larson’s relation on the column density, and 
attributed this to the fact that the clouds are in self-gravitational 
equilibrium. We plot CTy/r^/^ as a function of column density 
E = Mlnr^, where M is the mass and r is the radius of our 



Fig. 14: as a function of co lumn density a for the re¬ 

gions. In order to be consistent with iHever et al.l (l2009l) . CTy is 
the ID velocity dispersion as defined in Eq. |20]and E = Mlnr^ 
where M is the mass r is the radius. The solid straight line 
shows the boundary below which the structures are gravitation¬ 
ally bound. It is defined as CTy/r^/^ = . The gray 

cross es come from the catalogue of giant molecular clouds stud¬ 
ied in iRoman-Duval et H] (|201Q|) . 


regions (Fig. (H. Since we decompose one region into a set 
of nested subregions based on the G-virial, one single region 
will appear as a curve in Fig. M All our regions are quite close 
to being gravitationally bound, and the scatter around the y-axis 
is much smaller tha n what is found fo r Milky Way molecular 
clouds as studied in iHever et al.l (l2009l) . This difference arises 
mainly from the fact that the G-virial method tends to identify 
regions that are gravitationally coherent. 


7. Conclusion 

7.1. Results and perspectives 


The_£h^sicaLscaleJiere_canJiave_differen^nieanin2s^_s^H & 
Brun_^^ej^H200jijLazmaii_^^ogos^a|iL^ej£jj2004j^R”an-Duval 
et al. (^e.g. j2011h . 


In this paper we proposed a general method ( G-virial) for 
studying the structure molecular clouds. Different from the pre¬ 
vious methods which decompose molecular clouds based on 
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contours in the 2D plane or 3D PPV data cubes, in this method 
a map is generated by taking all the gravitational interactions 
between all the voxels in the 3D data cubes into account. 

The generated 3D G-virial maps have a dimensionless unit 
that is the same as that of the virial parameter, and a larger 
G-virial is related to a larger chance of being gravitationally 
bound. Therefore, the method provides a global picture of grav¬ 
ity in the PPV space. Using a hydrodynamical simulation model 
and a simple spherical-symmetric model, we demonstrate the 
connection between G-virial calculated in the observed PPV 
data cube and the G-virial calculated from the model. We 
found that the two are positively related. A larger G-virial 
is therefore linked to a larger chance of being gravitationally 
bound. Different from the virial parameter which quantifies self¬ 
gravity, the G-virial quantifies the effect of global gravity on 
the gas where all the interactions between all the particles in the 
3D data cubes have been taken into account. 

A map of this kind enables us to to identify regions in molec¬ 
ular clouds in terms of gravitationally coherent regions. Com¬ 
pared with previous ways of defining regions based on emission, 
the regions identified from the G-virial maps more coherent 
and are less dependent on the technical parameters. Another 
advantage of such a definition is that the regions identified are 
coherent under gravity and are likely to collapse on their own. 

We also demonstrated how to use the G-virial maps to 
quantify the structures of the identified gas condensations in 
terms of mass-radius and velocity dispersion-radius relations. 
We found that both the mass-radius relation and the velocity- 
dispersion relation of the model can be reconstructed with the 
help of the G-virial with a good fidelity when the model is 
symmetric. If the model is not symmetric, the G-virial map 
can used to quantify the structures of the gas condensations, and 
it allows us to compare different structures. 

As examples, we analyzed the ^^CO(l-O) emission from the 
Perseus and Ophiuchus molecular clouds, and found that both 
can be decomposed into several regions which are gravitation¬ 
ally coherent (Fig. n. Moreover, the cluster-bearing regions 
show higher values of G-virial at the centers which implies 
that gravity plays a more prominent role in these parts (Fig. [IB. 
We carried out an analysis of a total of five regions identified in 
the Perseus and Ophiuchus molecular cloud, and derived mass- 
size relations and velocity dispersion-size relations for the re¬ 
gions. We also found that cluster-bearing regions are more mas¬ 
sive at a given radius than those which do not bear a cluster (Fig. 
(IB, and the cluster-bearing regions have a higher velocity dis¬ 
persion at the outer parts (Fig. [TB . 

The method is general and can be applied to a variety of ob¬ 
jects observed in 3D PPV space where gravity is supposed to 
play a role, and from the observations the method offers ways to 
quantify and compare these structures. We leave these possibili¬ 
ties for further explorations. 

The importance of gravity in molecular clouds remains un¬ 
clear. The effect of gravity as compared to kinetic motion is usu¬ 
ally quantified using the virial parameter. The major uncertainty 
of the virial parameter comes from the definition of boundaries 
the regions. With the G-virial method, we can eliminate this 
uncertainty by focusing on the regions that are gravitationally 
coherent. As shown in Fig. M the G-virial-defined regions 
are much closer to being g ravitationally bound than the giant 
molecular clouds studied in iHever et (l2009l) . A study of a 
larger sample of molecular clouds can potentially tell if gravity 
is important throughout the Milky Way molecular clouds, and 
we leave this for a further study. 


7.2. Caveats 

Like all the other methods that quantify the ISM structures, our 
method is influenced by the line-of-sight confusion. One ad¬ 
vantage of our method is that the line-of-sight confusion can be 
reduced by taking the velocity information into account. For 
a single object, if the velocity dispersion increases with radius, 
the G-virial method can take this into account and reduce the 
line-of-sight confusion effect. If different objects have different 
velocities, they will be separated easily. If different objects with 
comparable mass velocity are found on the same line of sight, 
the G-virial method becomes inaccurate. 

The G-virial is a relative measure of the gravitational 
boundedness. As a result, it is not possible to tell whether a 
region is gravitationally bound based on the G-virial. How¬ 
ever, a larger G-virial value means a region is more likely to 
be gravitationally bound. It is also possible to derive the virial 
parameter afterwards based on the iso-G-virial contours. 
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Appendix A: Relation between G-virial and the 
virial parameter 

The relative importance between gravitational and kinetic en¬ 
ergy is usuall y characterized with the vi rial parameter, which is 
introduced in iBertoldi & McK^ (Il992l) . To evaluate the virial 
parameter, it is necessary to define a structure on which the virial 
parameter is calculated. 

One difference between G-virial and the virial parameter 
is that to evaluate the G-virial no such boundary is needed. 
Therefore the G-virial can be viewed as a generalization of 
the virial parameter to a continuous distribution of mass. 

To illustrate the physical meaning of the G-virial, we de¬ 
fine a quantity called particle virial, p-virial, which is the ratio 
between the gravitational energy of a particle to its kinetic en¬ 
ergy tTp-viriai = £^p/2£’k. The potential energy is determined 
by m X 0 where m is the mass of the particle and 0 is the grav¬ 
itation potential, and £’k is determined as 1/2 x m (v - Vc)^. To 
define the kinetic energy, the velocity of the center of mass Vc 
is needed. Therefore p-virial is only suitable in the cases where 
a center of mass can be easily found. In the case of a molecu¬ 
lar cloud, this is not straightforward because inside a molecular 
cloud the condensations can move at different velocities. The 
velocity of the center of mass depends on which clumps are in¬ 
cluded in the calculation, which is not unique. This difficulty is 
illustrated in Fig. lA.ll 
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Fig. A.l: The difficulty of finding a proper center of mass. We 
consider three clumps pi, p2, p3 moving at velocities vl, v2, 
v3. If we are interested in whether i is gravitationally bound, 
since the majority of attraction comes from pi, physically the 
center of mass should be the center of mass of pi instead of 
the center of mass of the three clumps pi, p2, p3, since p2 and 
p3 do not contribute much gravitational attraction to i compared 
to pi. Therefore to estimate the gravitational boundedness, a 
naive calculation of the center of mass where all the mass in 
the whole cloud is included is not appropriate. In the case of a 
complicated distribution of gas, finding the center of mass is not 
straightforward. 

One major advantage of the G-virial over the p-virial is 
that no center of mass is needed. The G-virial is the sum of 
the gravitational boundedness, and to evaluate it we do not need 
to introduce the center of mass. Therefore, the G-virial can 
provide estimates of the importance of gravity when the geome¬ 
try of the object is complicated and the center of mass is not well 
defined. 

Appendix B: Dependence on of G-virial the 
smoothing veiocity cq 

The only free parameter in our method is the smoothing velocity 
Co in Eq. [6l We added this parameter to avoid the divergence of 
Eq. [6] when the separation of velocity is zero. 

Physically, if the velocity separation is zero, and if two gas 
particles spatially coincide with each other, the gravitational in¬ 
teraction will cause them to collide. However, the velocity sepa¬ 
ration cannot be zero since the gas also has ve locity dispersi ons 
that have either thermal or non-thermal origin (iLarsonll 19811) . A 
physically lower limit of cq should be the sound speed. 

In reality, this is also affected by the superposition along the 
line of sight, since different gases that are physically unassoci¬ 
ated can stay along the same line of sight and therefore appear to 
be gravitationally bound. Therefore, it is difficult to find a unique 
value of Co based on first principles. In our case, we choose cq to 
be 1 km/s, which is larger than and still comparable to the sound 
speed. 

Here we investigate how our results can be affected by the 
parameter cq. First, we made a G-virial map by assuming 
Co = 1 km/s and then we made another map by assuming 
Co = 2 km/s. Finally, we made a voxel-by-voxel comparison of 
the maps (Fig. ED. We found the major effect of changing co 
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Fig. B.l: A comparison of G-virial map of the Ophiuchus 
molecular cloud under different cq. The horizontal axis is the 
G-virial calculated assuming co = 1 km/s and the vertical 
axis is the G-virial calculated assuming cq = 2 km/s. The 
grayscale image stands for the number of voxels that fall into 
each bin. 

on the absolute values of G-virial: if we change cq to 2 km/s 
the G-virial decreases by a factor of ~ 2. However, there is 
good correspondence between the old and new G-virial val¬ 
ues. Therefore, changing cq affects the absolute values of the 
G-virial; however, the relative values are unaffected. There¬ 
fore, all our figures should be unaffected by a change of cq and 
the G-virial axes of Fig. [T0]and Fig. [TTlwill change accord¬ 
ingly. 

In the case where the G-virial is applied to different data 
cubes, in order for the results to be comparable, a unique choice 
of Co is necessary. 

Appendix C: 3D renderings of G-virial of the 
Perseus and Ophiuchus moiecuiar ciouds 

In order to provide intuitive representations of our G-virial 
maps, w e present volume renderings of the 3D PPV data cubes 
with yt (iTurk et al.ll201 ll) in Fig. [71 A movie can be found in 
https: //vimeo. com/ll 5595991.. 

Appendix D: Comparison with other methods 

The previous methods to quantify the structure of molecular con¬ 
densations such as Clumpfind and Dendrogram focus on the 
structure traced by the intensity map. The Clumpfind method 
tends to produces isolated structures. The Dendrogram pro¬ 
duces hierarchical representations of the nested isosurfaces in 
3D molecular line data cubes, and offers multi-scale decomposi¬ 
tions. It can be applied to a much larger variety of situations. 
Both methods work in 2D (position-position) and 3D (PPV) 
space and both methods decompose the data to some extent. 

The G-virial method presented in this paper is also based 
on a map of intensity in 3D PPV space. However, its output is 
neither a list of structures (as in the case of the Clumpfind) nor 
a tree representation of a hierarchy of structures (as in the gas of 
the Dendrogram). Instead, it produces a map of the importance 
of gravity in the 3D PPV. The output is also a map 3D PPV space. 
The G-virial map itself does not provide a decomposition of 
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Fig. C.l: Volume rending representations of the ^^CO(l-O) 
emission and the corresponding G-virial map. The upper 
panel shows the result from the Perseus molecular cloud and 
the lower panel shows the result from the Ophiuchus molecular 
cloud. Blue stands for ^^CO(l-O) emission and red and orange 
stand for G-virial. Both clouds are projected along the veloc¬ 
ity direction. For the Perseus molecular cloud, the iso-surfaces 
of ^^CO(l-O) emission start from 0.3 K and increase in steps of 
0.64 K. The contours of the G-virial starts from 1.2 and in¬ 
crease in steps of 0.8. For the Ophiuchus molecular cloud, the 
contours of ^^CO(l-O) emission start from 0.3 K and increase in 
steps of 1.24 K. The contours of the G-virial starts from 1.2 
and increase in steps of 0.8. 
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the data. To decompose the cloud into regions and to analyze 
their properties, other methods are needed. 

In this work, we use the Dendrogram to identify gravita¬ 
tionally coherent regions from the G-virial map, and use the 
iso-G-virial contours to quantify the properties of the regions in 
the M-r and (Ty-r plane. The G-virial method is not a replace¬ 
ment of the other methods, but it provides a new map on which 
those methods could be applied. 
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